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1. Introduction 

In molecular dynamics simulations of heavy-ion collisions at intermediate energies 
(300 MeV/u - 2 GeV/u) it is still an open question question, which kind of 
forces should be used for the long range part of the nucleon-nucleon interaction 
which builds up the mean field. The Skyrme forces, which work very well for 
lower energies are not relativistically invariant. There are attempts to extend 
non-relativistic two-body potentials to higher velocities by requiring approximate 
Lorentz invariance [0, |^, ^, Q . Our aim is to derive from a relativistic lagrangian 
of the Walecka type a covariant Lagrange function for point-like particles which 
is suited for molecular dynamics and still has the saturation property of the 
original lagrangian. It should then be applicable in transport theories as for 
example QMD or BUU and should give reasonable results even at lower energies 
(50 - 200 MeV/u). 

In Section 2 the relativistic Hamilton function is derived in the small acceler- 
ation approximation from a many-body Lagrange function and it is expressed 
in terms of positions and canonical momenta of the particles. In Section 3 the 
QMD procedure is described, while in Section 4 we present our QMD calculations 
for fragment distributions in O + Br collisions at different energies and compare 
the calculated results with the experimental data. A few years ago the analogue 
calculation was done p with the Skyrme interaction, which is explicitly density 
dependent. It is interesting to compare those results with the new forces pre- 
sented in this paper in order to see how far they can describe saturation and low 
energy phenomena. 

The intention is, however, to go to intermediate energies, where QMD is more 
appropriate because quantum aspects like the Pauli principle or the uncertainty 
relation are expected to be less important. 

2. Relativistic equations of motion for nucleons 

In this section we are deriving the relativistic equations of motion for the mean 
4-positions and the mean 4-momenta of the nucleons. We start from a field 
theoretical lagrangian of the following type 



where we allow for two different ways to couple the scalar field If m*(0) = 
m — (7s0 we obtain the Walecka lagrangian, whereas for m*{(f)) = {m + g^cj))"^ 



+ 



5(vV^(x)7"V'(x)Aa(x) 
^A-{x){dpd^ + f,l)A^{x) , 
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we deal with the lagrangian proposed by Zimanyi and Moszkowski ^j. (In the 
following, if not expressed explicitly otherwise, we use units such that h = c = 1.) 
The field equations are 

{{td^ - g^A^ix)) - m* (0(x)) } = (2) 

for the nucleons, 

{df'dp + fxl) 0(x) = - ^m* (0(x)) ^(x)V^(x) (3) 

for the scalar field and 

(df'df, + ixl) = (?v^(x)7Xx) (4) 

for the vector field. 

In the mean-field approximation the source terms are replaced by their expec- 
tation values ip{x)4'{x) =^ Ps{x) and '?/^(x)7°''?/'(x) =^ j'^ix) so that the scalar and 
vector fields become classical fields. For classical fields, however, an arbitrarily 
weak time dependence in the scalar density Ps{x) or current density j°'{x) results 
in a radiation field which travels away from the nucleons. The total field energy 
in the radiation may even be less than the masses of the field quanta. Since the 
fields are only effective fields, and for the scalar field there is no corresponding 
elementary particle, the physical meaning of this radiation is obscure. 

As we are using the relativistic scalar and vector fields only for the mean-field 
part of the nuclear interaction we follow the suggestion [|], || to exclude radiation 
of fictitious mesons from the very beginning by using the act ion- at- a- distance 
formulation with the symmetric Green's function of Wheeler and Feynman . 

Gs,v(x -y) = l {Gt'~\x -y) + - y)) (5) 

The formal solution of the wave Eq. (|^) (for m* = m — gs(f)) 

(Pix) = Qs i d^y Gsix - y)i'{y)ij{y) (6) 



andofEq. (D 

A-{x) = I d% G^{x - 2/)^(y)7"^(2/) (7) 



fulfills the desired boundary condition of vanishing incoming and outgoing free 
fields. 

Eliminating the fields from the lagrangian (^ leads to the non-local action 
which contains only nucleon variables: 



d^x C{x) = j d'^x ^(x)(7"z9« - M)ij{x) 

+ l^al I d'^x dS i){x)i){x) G^{x - y) i){y)i){y) 



]^gl / d^x d^y V'(x)7"V^(x) G^{x - y) il){y)-ic,'^{y) 
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If quantum effects are negligible the scalar density and vector current density can 
be represented in terms of the world lines of the nucleons as 



P3(x) := ^(xMx) E /rfTj 5^(x-rj(rj)) (9) 

r{x) := ij{x)ri:{x) ^ j^fdr. 6\x - r^{r. )) u^{t. ) , (10) 

where rj(rj ) denotes the world line of the nucleon j at its proper time Tj and 
«j(rj ) its 4- velocity. Nucleons in nuclei can however not be localized well enough 
for treating them as classical particles on world lines. The Fermi momentum sets 
a limit to the radius of the nucleon wave packet in coordinate space which is of 
the order of 2 fm. These effects, the Pauli principle and other quantum effects, 
are properly taken care of in Antisymmetrized Molecular Dynamics (AMD) ||Tl| 
and in Fermionic Molecular Dynamics (FMD) [|12]. However, presently neither 
the AMD nor the FMD equations of motion can be solved numerically for large 
systems like gold on gold. Therefore, we mimic the finite size effect of the wave 
packets by folding the finally obtained forces with gaussian density distributions 
for the nucleons. This "smearing" provides also the prescription how to eliminate 
strong forces at short distances which are taken care of by the random forces of 
the collision term. 

Therefore, the forces to be derived in the following are only meant to represent 
the long range part of the nuclear interactions which constitute a mean field. In 
this sense rj{T- ) is the mean world line of the nucleon and Uj{T- ) its mean 4- 
velocity. 

Using the representation (^ the scalar field at a space-time point x, for ex- 
ample is given by integrals over past and future proper times Tj as 

m=9sj: I dr. G,{x-r;{T.)) . (11) 

By identifying the following expressions with their representation in terms of 
mean positions and mean momenta of the nucleons 



d'^x d'^y ^{x)ijj{x) Gsix - y) ^{y)tp{y) 

Y: Idndr^ G^{r,{r,)-r,{r^)) , (12) 
i,j=i 

d^x dS tp{x)'j°'i){x) G^{x - y) tp{y)'jai'{y) ^ 

X: / dr, dr. GMir, ) - r,{T^ )) uUr, )uf{T^ ) (13) 



i,j=l 
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in the action one obtains the non- instantaneous action-at-a-distance 
d'^x C(x) 



i=i •' 

+ \glj: jdr.dr^ G,{n{r,)-r,{r^)) (14) 

i.j=l 

- 2^^^ J G'v(ri(ri ) - rj(rj )) u^^ir^ )uf (tj ) . 

i.j=l 

Lagrange multiphers Ai(rj ) are introduced to ensure Mi(rj )^ = 1. 

The equations of motion which result from this Wheeler-Feynman action are 
equivalent to the classical field equations, provided the system does not radiate 
0. 

However, the Wheeler Feynman equations of motion cannot be solved in gen- 
eral as one needs to know the world lines for all the past and future times in 
order to calculate the fields which enter the equations for the world lines 
Furthermore, the no-interaction theorem states that, except in the case of non- 
interacting particles, there exist no covariant equations of motion for world lines 
in which only the 4-positions and the 4-velocities at a given time enter, as is the 
case in non-relativistic mechanics. 

In the following this no-interaction theorem is circumvented by introducing an 
approximative solution to the non-local Wheeler Feynman equations of motion. 
This is achieved by the so called small acceleration approximation which does 
not assume small velocities. 



2.1. Action-at-a-distance in the small acceleration approximation 

In order to conserve manifest covariance of the equations of motion we intro- 
duce first the concept of a scalar time. For that the whole Minkowski space is 
chronologically ordered by a set of space like surfaces S{x) which attribute to 
each 4-position x a scalar time t by 

S{x) = t , with d^S{x) > , d°d°S{x) > . (15) 

The simplest choice for these isochrones, which we shall take in the following, is 
S{x) = rjax"", where rja is a position independent time like 4-vector. 

Unlike in a collision the nucleons are not strongly accelerated by the action 
of the mean-field. Therefore, in order to describe the motion in the mean-field 
one may expand each world line around the proper time t-^ at which the particle 
is at a given isochrone, i.e. t = S{rj{Tp) for each j. 
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rfir, )=rfirp + {r^- + ^(t. - a^ir^ + ■ • • . (16) 

This way we are defining a synchronization prescription for all particles, which 
does not depend on the frame. For calculating the fields we neglect of the 
quadratic term with the acceleration (if{Tp and all higher powers which leaves us 
with a straight world line in the vicinity of the synchronizing time t = S{rj{Tp) 
(c.f. Fig. 1). This is called "small acceleration approximation" 

The small acceleration approximation is of course best in the vicinity of r^{Tp 
and becomes worse further away. As sketched in Fig. 1 only those parts of 
the world lines (thick lines), which are outside the light cone (centered at x), 
contribute to the field strength at point x. A world line which hits the light cone 
far away from x may be badly approximated by Eq. (|16|) , but for short range 



interactions a distant particle does not contribute anymore to the field at x. 

Thus, the first condition for the validity of the approximation is that the range 
fi~^ is small compared to the curvature of the world lines, i.e. the inverse of the 
acceleration. The second condition is weak radiation, which is fulfilled when the 
acceleration is small compared to the meson mass fi. Both conditions are actually 
the same, namely 

I fl/.a'' I < /i^- (17) 

The assumption of small accelerations is justified if the (p and A" fields are only 
meant to be the mean-field part of the nucleon-nucleon interaction in a hadronic 
surrounding. The hard collisions between individual nucleons which are due 
to the repulsive core will cause large accelerations and also create new particles. 
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These hard colhsions cannot be described within Walecka type mean-field models. 
Therefore, it is consistent to regard and as Hartree mean- fields which 

bring about only small accelerations and which are not radiated away from their 
sources. 

Inserting the straight world line into Eq. (^) results in the easily understood 
situation that the field at a point x is just the sum of Lorentz-boosted Yukawa 
potentials which are traveling along with the charges: 



where Rj{x)'^ is given by 



47r Y 



exp |-/is-y/-Rj(x)2| 



R,ixf = -ix-r,{TpY+ ( 



X 



The vector field is derived in an analogue fashion as 



X 



(19) 



(20) 



At this level of the approximation the causality problem with the advanced part 
of the Green function is not present because the retarded and the advanced fields 
are identically the same when they are created by charges moving on straight 
world lines. Therefore one may regard the fields as retarded only. In addition 
the unsolved problem of radiation reaction, where the radiation acts back on the 
world lines [113], does not occur because there are no radiation fields any more. 



2.2. Instantaneous action-at-a-distance 



In the spirit of the small acceleration assumption discussed in the previous sec- 
tion one can use the straight line expansion in the action (|1^ and perform the 
integration over Tj . This results in an instantaneous action-at-a-distance where 
the Lorentz-boosted Yukawa fields appear again and there is only one time, the 
scalar synchronizing time t. 



- d^S{n{t))ur{t) 



-{M-Ut))u,{tf 



^ l^gle^v{-^^s^Ri^{tf] 



2 ^47r 
J 



(21) 



_l^exp{-^v£^ 
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The four-positions rp(t) = rf^r^ (t)) and 4- velocities uf{t) = uf{r-^ (t)) are to be 
taken at the same scalar synchronizing time t and 

^SW := -{rm - r,{t)f + - r,,(t))u^{t)\" (22) 

The small acceleration approximation together with the introduction of a syn- 
chronizing hj^ersurface S{x) leads to an equal time lagrangian which is Lorentz- 
scalar and written in a manifestly covariant way. 

Giving up explicit covariance and choosing rj = (1,0, 0,0) in a special coordi- 
nate frame, the positions and velocities take the form 

r,{r, it)) = (t,m) and u,{t, [t)) = ^ (l , m) ■ (23) 

With that a Lagrange function C{ri{t) , Vi{t)) can be defined which depends only 
on the independent variables and one time. Even the Lagrange multipliers Ai(t) 
are not needed anymore if the variation is with respect to V[{t) instead of all four 
<(t). 

The advantage of the instantaneous Lagrangian is that one can define easily 
the hamiltonian and the total momentum, which are then strictly conserved by 
the equations of motion. 



2.3. Hamilton equations of motion 

In the following we want to express the total hamiltonian as a function of the 
positions ri and the canonical momenta pi. Using relations (pTD-(p3D the Lagrange 
function is 



E 



m 1 



lib i X ^ 

^^2^ Anu? 



gl exp{-/Xs/R| 



J 



1^ exp {-^^J Rl{t)\ 



2^ 4™? 
J ' 



(24) 



The canonical momenta and the energy are determined from the Lagrange func- 
tion above as 



Pi 
E 



dC 

dvi 



(25) 
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Performing these derivatives, the energy is given by 
E 



2 W 



0\2 



[U, 



0n2 







(26) 



+ 2 ^ ^^^j^'^ 
where we use the abbreviations 



0n2 



fl-v exp{-/ivRij} 



RjidR/'~~'" 4:71 
Rl = r-^ + (fyWj)^ 



exp{-/isRji} 



'J ' 



(27) 
(28) 



and ^ji is defined in the same way as /ji. It is worthwhile to note that /y 7^ f]i. 
The momentum of a particle looks as 



Pi 



mui — 



1 



E k 



1 ^ 

- 9 E /jiT^ [^ij(^ij^i) + ^i(^ij^i)' 
J J 



+ ^ E ( + z^9]i ) - ^ E 



2 ^ 



+ 2^ 
j 



U 



2 ^ 

J 



Mi 



(29) 



rij(fijMi) +Mi(fijMi) 



Let us first prove that for isotropic nuclear matter, like in the original mean- 
field picture [jl5| the vector potential will not contribute to the canonical mo- 
mentum. In that case the summation can be replaced by an integration over 
space and momentum, folded with the phase space distribution function. Since 
for isotropic nuclear matter the distribution function is independent of the posi- 
tion, the space integral can be carried out immediately and the remaining part 
is written as a summation over momenta. This way we get for the Yukawa terms 
in nuclear matter 
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E h ir,u-y ^ y: I d'r, KR^{f^;f = (l) ^ E 1- (30) 



Substituting the above expressions into the momentum (|2^) the /ij term gives 
the same contribution as the /y, while the terms cancel each other and the 
momentum can be written as 

The last term gives zero, since the average of mj is zero in nuclear matter and 
Eq. (^) can be written as 

Pi = m*ui , (32) 

where 



171=171 




u 



I ■ ■ I ^ /j72*2 _|_ 



Similarly, the energy density (|26|) has the form 

/ \ 2 



e = E V"^*^ + + ^ (^^ j + ^ j ("^ - ■ (34) 

Eqs. ( P3| - P^ are just the mean-field results of the Walecka model [^]. 

It is worthwhile to mention that the small acceleration approximation which 
introduces instead of | | in Eqs. ( p7H28| ), is definitely needed to get back 



the relativistic mean-field result for nuclear matter. Therefore, for the coupling 
strengths Qs and we can simply take the values obtained from the saturation 
properties of the nuclear matter [^]. 

For the Hamilton equations we need the energy as the function of the canon- 
ical momenta instead of the 4- velocities u\. For that we have to invert the p{u) 
equation (p9D into a u{p) one. One can do this approximatively by first expand- 
ing the 4-velocities in the p{u) relation up to the third order in the quantities 
{p/m)'^, /ij and ^fy, and then substituting that expression for u{p) into the energy 
Eq. (P^l). At the end we will check the validity of this approximation. After 
tedious calculations the energy can be expressed as follows (in the following we 
include explicitly c and Ti) 
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E 



2E ^ji 

iJ 



+ .E 



E f0i - E ^ji^j 



(35) 



where 



n \ exp{-^^ij} 



mc 



and 



~ ^2 



(36) 



with pj = pi/mc. Qij is defined similarly to /ij. The last term is already of third 
order and hence small. 

From Eq. (^) the Hamilton equations 



d ^ BE d ^ dE 



dt 



r; = - 



(37) 



are calculated easily. 



3. Details of the QMD calculations 



Molecular dynamics is a classical many-body theory in which some quantum 
features due to the fermionic nature of nucleons are simulated. For the details 
of the theory we would like to refer the reader to the works of Aichelin [|16| and 
the Frankfurt group |T^. In our calculation we introduced some modifications 
which are described in |l^ . In the following this model is used in a relativistic 
treatment with scalar and vector forces. 



3.1. The mean-field forces 

The scalar and vector forces are the relativistic two-body forces derived in the 
previous section which give saturation for nuclear matter even without an explicit 
density dependent term. However, Walecka's values for the coupling constants 
result in a too high compressibility. For this reason it seems to be more adequate 
to use instead the Zimanyi-Moszkowski (ZM) lagrangian P] which provides a more 
reasonable compressibility. The different coupling of the mean-field (m* = (m + 
gs4')~^) modifies the scalar part of the total energy. In the used approximation, 
the expression 
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11.04 


13.74 


7.0 


550 MeV 


783 MeV 


ZM 


7.84 


6.67 


7.0 


550 MeV 


783 MeV 



Table 1: Parameters of the Walecka and the Zimanyi-Moszkowski forces 



2 ^ 



h 



Pi 



1 + 2E /j 



V 



k 



jk 



(3^ 



replaces the second term on the r.h.s. of Eq. (pS]). Here we see the effect of the 
nonlinear coupling in the ZM lagrangian. The strength of the scalar potential 
which is felt by particle j is weakened. 

We compare the Walecka forces with the ZM ones and study their effect on 
nuclear multifragmentation. The coupling strengths Qs and are determined 
from the ground state properties of nuclear matter [0, [T^. In order to include 
the symmetry energy we take also the g meson into account. To simplify the 
calculations the mass of the g meson is taken to be equal to the one of the u 
meson, and the coupling constants are fitted to get the appropriate symmetry 
energy. In Table 1 we list all coupling constants used. 

Since the relativistic two-body forces are meant to describe the long range 
part of the interaction, the Yukawa functions are folded with a gaussian density 
distributions for the nucleons. For nuclear matter one gets a smooth density 
distribution if the gaussians e""" ('"-'^i) have a width parameter of a = 0.5 fm~^. 
We use that value for folding the Yukawa forces, analogue to our earlier QMD 
calculations |l8l, with Rf^ replacing r^^ everywhere. 



3.2. Initial conditions 

For the initial coordinate and momentum distribution of the particles we use the 
ones, which were prepared for the non-relativistic calculations P . The main 
aspect when generating initial positions and momenta is to get a smooth phase- 
space distribution. The ground-state energies calculated with distributions which 
fit nuclear charge densities are found to be good for both, the Walecka and the 
ZM forces. 
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Figure 2: Distribution of nucleon-nucleon relative momenta up to time t in 
Ca+Ca central collision 



3.3. Cross sections 

The multifragmentation depends on the cross sections used. According to Cugnon, 
his parameterization [|l^ gives a good fit only for relative momenta > 0.8 GeV/c 
in the two particle center of mass frame. For Ca + Ca central collisions at 200, 
400, 600 and 800 MeV/u initial beam energies Fig. 2 shows however, that the 
vast majority of the collisions occur at relative momenta below 0.2 GeV/c. To see 
the effect of the cross section on multifragmentation we fit the experimental free 
nucleon-nucleon cross sections as a function of the relative momentum with 
second order polynomials and use that fit in our calculations. However, due to 
the screening effect of the other nucleons we do not allow collisions at distances 
larger than Tcou- The screening length is somewhat arbitrary, so we work with two 
parameterization, rcou = 1.6 fm and rcou = 2.4 fm. The low energy cross section, 
we use in the calculation, is displayed in Fig. 3. The dependence of the fragment 
distribution on the different cross sections is discussed in the next section. 



3.4. Effective mass in two-body collisions 

In relativistic many-body calculations single-particle energies Ci cannot be defined 
such that (ei,pi) form a 4-vector. This holds true only for the total momentum 
P and total energy E. However, we still try to deduce a single particle energy, 
which we use for prescribing the conservation of energy in each collision. If one 
considers only the scalar potential in the lowest order the solution is easy. The 
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Figure 3: Low energy nucleon-nucleon cross section 



equation (|35D in the leading order can be written as 
E 



^ +^E^jiVl+^i (39) 



2tr i+2E/j.yr7f 2 

There are two ways in our case to define the effective mass m* appearing in 
the coUisions: 

Case A: If we neglect the momentum dependence in /y and g-^, the energy of 
the i*^ particle can be written as 

E{A)-E,{A-l) 



1,-2 1 ^ /ij 

1+Pi-oE 



1 



2t ^/^Tf i + 2E/.k 



(40) 



1 ^ / ^2 ^ 

+ 2 E 5'jiV 1 +Pi + terms independent of p; , 
j 

where Ei{A — 1) means the energy of the system when particle i is taken out. 

-,2 

We can expand Eq. P0| ) in and get 
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Figure 4: The evolution of the fragment distribution in time in 200 MeV/u O + 
Br central collisions 



where eio is independent of p^. This can be written as 



where the effective mass mf turns out to be 



1 



771; 



/ij 



1 



2r l+2E/.>^ 2 ^ 



(41) 



Case B: If we assume, that the momentum dependence of and g^i is the 
same as in nuclear matter, that is 



fin. 



J^i ^ V J k JT^i 



+ - J: gir,) , (42) 



[see Eq. (^)),then the effective mass for finite systems turns out to be 

/(^ij) 
^ l + 2^/(%) 



rri: 



1-E 

j 



(43) 



k 

kT^j 



We found in our calculations that the definition (^T]) gave better conservation of 
the total energy. 
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Figure 5: The evolution of the fragment distribution in time in 50 MeV/u O + 
Br central collision 

3.5. The treatment of a relativistic collision 

The individual two-body collision is calculated in a relativistically covariant way, 
as given in 



21[] . The particles are allowed to collide if the impact parameter 
b = \l Ri2 — hi2/vi2 with 

2 , fMx2-X^)Y 



(44) 



R 



12 



h 



12 



-(Xi - X2) + 

Y rriic ) 

Pl{x2-Xi) P2{X2 - xi)mic 



rriiC 

I mim2C^ 
V P1P2 



P1P2 



and 



is smaller than the collision radius ycr/vr. In the model the collision takes place 
at equal time in the CM of the total system, provided the closest approach of 
the world-lines occurs within the following time step. After the collision we only 
allow final momenta which do not let the particles approach each other further. 

3.6. Pauli blocking 

Pauli blocking is treated as in Ref. [Q. We calculate the averaged phase-space 
distribution function at the location of each particle. The collision can occur 
if the new phase-space distribution function is smaller than 1 for the colliding 
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Table 2: The contribution of the different energy terms of Eq. ( pS] ) in MeV per 
particle for gold on gold systems for very extreme cases. The upper numbers 
refers to the Walecka, the lower ones to the ZM forces. The calculations are 
made for 1 MeV/u and 2 GeV/u energies at normal and 3 times normal nuclear 
matter density. 

particles and for all other particles after the collision. Since the distribution 
function at the phase space location of a particle depends on the momenta of 
all the other particles, without this condition it may occur that the phase space 
density becomes greater than one even for particles which did not take part in 
the collision. 

3. 7. Definition of fragments 

We consider the particles to belong to the same fragment if they stay together in 
coordinate and momentum space and are bound. In our calculation we follow the 
evolution of the system up to 360 fm/c for higher, and up to 720 fm/c for lower 
energies. We see in the Figs. 4 and 5 that in the last 120 fm/c the so determined 
fragment distribution does not change anymore. 

In order to calculate the energy of a cluster we calculate the /y and gij in the 
total CM frame because they are Lorentz-scalars. Only the momenta appearing 
in the momentum dependent terms have to be Lorentz transformed to the center 
of momentum frame of the cluster. The radii and the energies of the fragments 
are very near to the experimental values, although some fragments can shrink a 
little. However, we do not get energies lower than —8.8 MeV/u for any case. 

4. Multifragmentation in the O + Br collision 

Among the experiments carried out to study nuclear multifragmentation, emul- 
sion measurements were the first to give an almost exclusive identification of the 
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0+Br central collisions (500) events/energy 





5 10 15 20 25 30 

Z 



Figure 6: Fragment distribution in O + Br central collisions at energies 50, 75, 
100, 150 and 200 MeV/u 



atomic numbers of the fragments emitted in heavy ion collisions. The investiga- 
tion of the same system at several bombarding energies is particularly useful for 
testing models. Such data are available for the O + Br system at bombarding 
energies 50, 75, 100, 150 and 200 MeV/u [^. Nowadays modern experimental 
techniques allow mass production of such measurements ||23|. Few years ago we 
carried out calculations for the emulsion experiment with a non-relativistic force 
to test the treatment of the Pauli blocking ||^. Now we have made the analogue 
calculation with relativistic forces at different energies to check their applicability 
at lower energies. These forces are meant to be used at higher energies, but they 
should also have good non-relativistic properties. 

As a first step we checked the validity of the expansion in /y and g^. In 
Table 2 we give the values of different energy terms of Eq. (|3^) for a gold on gold 
system which was artificially compressed to different densities and given different 
relative energies. One can see, that for the ZM case the contribution of the meson 
fields are much smaller than the nuclear rest mass mc^. Even for the Walecka 
case, except maybe for the high compression with j° = 3go, the expansion is 
applicable. Thus our approximation in the expansion (|29|) of the momentum is 
justified. 

In Fig. 6 the fragment distributions of the O + Br collisions are displayed for 
different energies. The contribution from the largest fragments are very small. 
Due to the computer time needed, the figures present results averaged over only 
500 events. Thus the statistical fluctuation for events with small yields are still 
large. In Fig. 7 we compare the results obtained with the Walecka and ZM 
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O+Br 200 MeV/u 
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Figure 7: Fragment distribution in O 
different cross sections and forces 



Br central collisions at 200 MeV/u for 



forces using the Cugnon parameterization [|l^ and our fitted cross sections for 50 
MeV/u and 200 MeV/u bombarding energies. The ZM force gives a somewhat 
better agreement with data than the forces derived from the Walecka lagrangian. 
The effect is larger for higher energies. However, using the Cugnon cross section 
gives more or less the same results as using the one fitted to the low energy free 
nucleon-nucleon scattering. 

The fragment distributions shown in Fig. 6 do not differ significantly from the 
earlier ones calculated with non-relativistic Skyrme forces, however one should 
mention, that for lower energies the experimental cut does not correspond to the 
central events. In Fig. 8 we display the distribution of impact parameters for 
the 10% highest multiplicities. For the 200 MeV/u collision the multiplicity cut 
selects the central events (impact parameter is less than 3 fm), but for 50 MeV 
it results in broad distribution. The yields for 50 MeV/u and 200 MeV/u beam 
energies with the multiplicity cut are given in Fig. 9. 



5. Conclusion 

From an effective field theory lagrangian with scalar and vector mesons we de- 
duced in the small acceleration approximation coordinate and momentum de- 
pendent relativistic forces, which are suitable for transport equation calculations. 
Even without an explicit density dependence, saturation of nuclear matter is ob- 
tained at proper densities and energies due to relativistic effects. The Zimanyi- 
Moszkowski forces, where a derivative coupling is used, can produce saturation 
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Figure 8: The distribution of impact parameters to the highest multiphcities. 
The sohd hne corresponds to 50 MeV/u colhsions and selects 11.3% of the events, 
while the dashed line represents 9.5% for 200 MeV/u 



even without relativistic effects. 

The proposed forces will be used in the intermediate energy domain (500 
MeV/u - 2 GeV/u), where the mean- field plays still an important role and rel- 
ativistic forces are necessary. However, they should work at lower energies too. 
As a first check we studied their low energy behavior in a QMD calculation. We 
investigated both, the more static features of these forces in the ground-states of 
nuclei and their dynamical properties in the fragmentation process in heavy ion 
collisions. The fragment distributions in the O + Br (50 - 200 A MeV) collision 
were calculated. The result we obtained is of similar quality as the one with the 
non-relativistic Skyrme forces M, ITH . 



That encourages us to use these forces for beam energies in the domain of 
500 MeV/u - 2 GeV/u, where non-relativistic Skyrme type forces cannot be used 
anymore. A further test of the model should be for example the comparison of 
the calculated flow (which is especially sensitive to the momentum dependence) 
for different energies with experimental results of the FOPI collaboration ||23|| . 
Such calculations are under progress. 
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Figure 9: The yields for 50 MeV/u (upper) and 200 MeV/u (lower) beam energies 
with trigger to the events with the 10% highest multiplicities. The solid line 
displays the yield obtained for a central collision (b=0.5 fm), the dashed line 
displays the average over the highest multiplicity events. The experimental data 
are shown by crosses 
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